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Bose-Einstein condensates of dilute gases are well-suited for investigations of vortex dynamics 
and turbulence in quantum fluids, yet there has been little experimental research into the approaches 
that may be most promising for generating states of two-dimensional turbulence in these systems. 
Here we give an overview of techniques for generating the large and disordered vortex distributions 
associated with two-dimensional quantum turbulence. We focus on describing methods explored 
in our Bose-Einstein condensation laboratory, and discuss the suitability of these methods for 
studying various aspects of two-dimensional quantum turbulence. We also summarize some of the 
open questions regarding our own understanding of these mechanisms of two-dimensional quantum 
turbulence generation in condensates. We find that while these disordered distributions of vortices 
can be generated by a variety of techniques, further investigation is needed to identify methods for 
obtaining quasi- steady- state quantum turbulence in condensates. 



1. Introduction 

Dilute-gas Bose-Einstein condensates (BECs) present unique opportunities for studying quantum 
fluid dynamics due to the wide range of experimental tools available for probing and manipulating 
condensates, and the quantitative accuracy available with theoretical and numerical approaches 
[l]. Experimentally, magnetic and laser fields can be arranged to create a wide variety of trapping 
configurations and geometries. Superfluidity hallmarks such as persistent currents and quantized 
vortices [2| [5] may be precisely generated and their dynamics probed, and numerous measurement 
techniques allow for observations of density and quantum phase distributions. Regarding simulations 
of superfluid dynamics in BECs, common numerical approaches utilize the Gross-Pitaevskii equation 
(GPE), a nonlinear Schrodinger equation that incorporates a mean- field approximation to represent 
the BEC with an order parameter, the familiar macroscopic condensate wavefunction 4 . Numerous 
variants of the GPE such as those that model finite-temperature environments extend the predictive 
accuracy of simulations [5J . Taken together, these features allow for highly synergetic experimental, 
numerical, and theoretical explorations of superfluid dynamics. 

In this respect, Bose-Einstein condensates may be ideally suited for investigating quantum turbu- 
lence and the dynamics of interacting quantized vortices in compressible quantum fluids; for recent 
reviews of this subject see Refs. 6 and [7 and the references therein. In superfluids, vorticity ap- 
pears in the form of superfluid- free vortex cores about which there is quantized flow circulation [2]. 
Whether in superfluid helium or trapped BECs, quantum turbulence is often phenomenologically 
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described as a tangled distribution of these quantized vortices throughout the fluid [8] [9] , and recent 



work with superfluid helium provides fascinating visual evidence for these vortex distributions 10 
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In BECs, quantum turbulence is relevant to early investigations on the formation of large vor- 
tex lattices [15] and routes to vortex nucleation by dynamical instabilities 13] fl4l, and has been 
experimentally observed and explicitly addressed in these contexts 
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More recent 

experiments have been principally dedicated to the study of quantum turbulence phenomena and 
the related vortex distributions [l9j [20J [2l] [22] . Images of vortex tangles in BECs have also been 
experimentally obtained [l7|[l9]. Yet while BECs readily permit vortex observations, other aspects 
of turbulence remain significant experimental challenges within the BEC field. Notably, direct mea- 
surement of kinetic energy spectra is an open problem, whereas such measurements of liquid helium 
have allowed direct comparison of classical and quantum fluid turbulent energy spectra [23 . The dif- 
ferent approaches and systems for studying quantum turbulence are therefore complementary, with 
BEC research playing a potentially important role in understanding the broad subject of turbulence, 
particularly in the context of compressible quantum fluids. 

More specifically, BECs are now beginning to contribute to the development of an under- 
standing of two-dimensional (2D) turbulence 
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in compressible quantum fluids. In the 
basic manifestation of 2D classical turbulence (2DCT), energy and enstroph^] are injected into 
a fluid at a small length scale, and patches of vorticity merge and lead to the development of 
large-length-scale flow. This is opposite to the dynamics of the familiar energy cascade of three- 
dimensional (3D) turbulence in a classical fluid, in which large-scale flows continuously decay to 
smaller-scale flows that are eventually dissipated by viscous damping. This inverse energy cas- 
cade of 2DCT has been studied for decades and numerous excellent reviews exist on 2DCT 27 28 
29 . 



1.1. Experimental study of two-dimensional quantum turbulence in BECs 

In order to study two-dimensional quantum turbulence (2DQT) in a BEC, one might envi- 
sion subjecting the BEC to continuous injection of energy and vorticity. Dissipation of en- 
ergy may also exist, perhaps by thermal damping that could induce vortices to exit the BEC 
at its boundary, or perhaps due to vortex-antivortex annihilation. A balance between forc- 
ing and dissipation can be envisioned to lead to a statistically steady degree of turbulence in 
a BEC. We might picture such a state as an approximately constant mean number of vortices 
in the BEC, since enstrophy in a quantum fluid is proportional to the number of vortices 
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With the likelihood of future experimental methods for observing vortex dynamics and mea- 
suring energy spectra, we further imagine a scenario in which 2DQT can be experimentally and 
numerically studied, linked with new analytical approaches, and compared with 2DCT phenomena. 
The prospects for the development of a new understanding of at least this one aspect of the excep- 
tionally challenging topic of turbulence are tantalizing motivations for pursuit of 2DQT research. 

Unfortunately there are a number of experimental challenges that must first be overcome. Sig- 
nificantly, an isolated, trapped BEC subjected to continuous injection of energy would eventually 
become depleted of atoms as the system rethermalized. It appears that a BEC experimentalist's 
primary hope in this respect is to develop methods for studying 2DQT that are quasi-steady, where 



a Over a 2D surface S, the enstrophy Q = J s dS \ou\ 2 is a measure of the vorticity uo = [V x v] z about the z axis, 
where v is the fluid velocity in the plane normal to z. 
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steady forcing and dissipation rates are balanced over finite timescales during which the BEC is not 
significantly depleted of atoms. Experimentally, the path to studying decaying 2DQT is perhaps 
much simpler, since this topic is more concerned with the dynamics of a BEC after forcing has 
stopped. 

While it should not be surprising that disordered vortex distributions in a highly oblate BEC are 
relatively simple to obtain by a variety of experimental techniques, as we discuss in this article, how 
best to obtain a quasi-steady state of forced 2DQT remains an open problem. We have begun to 
empirically tackle only the forcing aspect of this problem, although simulations are also beginning to 
address this issue 22 31 32 33 . Several of our methods for driving large numbers of vortices into 



highly oblate BECs are described below. We generally associate these states with 2DQT in the sense 
that the disordered vortex distributions may be considered a 2D phenomenological equivalent of the 
vortex tangle that exists in 3D quantum turbulence. We do not provide any evidence that such dis- 
ordered vortex distributions conform to particular kinetic energy spectra or display vortex dynamics 
that may be expected in fully developed 2DQT; discussions of links between vortex distributions, vor- 
tex dynamics and aggregates, and kinetic energy spectra in BEC 2DQT can be found elsewhere 
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31; 33 34 35 . We also note that in the highly oblate trapping limit, BEC confinement in one 
spatial direction is much stronger than the trapping in the two orthogonal directions. This limit is 
not strictly necessary for studying 2DQT, but vortex excitations in the form of Kelvin waves are 
suppressed and vortices can be accurately approximated as having point-like dynamics in a plane 
within this limit [36] . 

In the following, we first describe some of the key phenomenological features of 2DCT, and 
summarize only a few of the theoretical advances in 2DQT, some of which have been aimed at 
uncovering the similarities and differences between the classical and quantum systems. This article is 
not intended as a review of the state of theoretical understanding of 2DQT in BECs; for discussions 

and the 
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of theoretical aspects of 2DQT we direct the reader to Refs. 
references contained therein. The study of 2DQT, particularly in compressible quantum fluids such 
as BECs, is still young and a full range of characteristics remains to be explored. Our aim is to give 
an overview of the types of observations that may be desirable in 2DQT BEC experiments, then 
to describe our experimental approach to creating and probing 2DQT in highly oblate BECs. We 
present a range of tools available for exciting the BEC using laser beams and magnetic fields. The 
bulk of the remainder of this article is devoted to describing various specific methods for driving 
vortices into a highly oblate BEC, and we show example images of BECs in such excited states. 
We conclude by briefly summarizing some of the open experimental challenges in generating and 
probing 2DQT in BECs. 



2. Overview of Two- Dimensional Turbulence 

From a phenomenological standpoint, Lesieur describes classical turbulence of continuum flows as 
having non-deterministic flow details, rapid mixing relative to molecular diffusion, and flow charac- 
teristics distributed among and interacting over a wide range of length scales [38]; see also Sommeria's 
similar characterization of turbulence 27 . Such a description encompasses 2D turbulence, where 
flow dynamics vary primarily over two spatial coordinates [x and y); flow in a third direction (z) 
may exist uniformly throughout the x-y plane or there may be negligible flow along z, hence a 3D 
fluid system may display 2D turbulence. Building from this general description, the further char- 
acteristics of 2D turbulence are markedly different from the characteristics of 3D turbulence. As 
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noted earlier, one such difference is the direction of energy flux between length scales: in 2D flows, 
small-scale forcing can lead to the development of large-scale flows. 



2.1. Key ideas of 2D turbulence 

Given their disordered and unpredictable nature, and challenges in discerning details of flow dynam- 
ics, turbulent flows are often characterized by the statistical properties of the system such as the 
kinetic energy density E(k) over wavenumber /c, which describes the distribution of kinetic energy 
among the length scales of the system. In 1967, Kraichnan 24 found the existence of two inertial 
ranges associated with distinct kinetic energy spectra power laws, joined by the wavenumber asso- 
ciated with energy and enstrophy forcing. The first range extends from the forcing length scale to 
larger length scales (or smaller wavenumbers), and displays a kinetic energy spectrum E(k) oc k~ 5 ^ 3 . 



This spectral distribution has the same form discovered by Kolmogorov in 1941 39 in the context of 
3D fluids, but in 2D this spectrum corresponds to an inverse energy cascade. The second range ex- 
tends from the forcing scale to smaller length scales (or larger wavenumbers) and the kinetic energy 
spectrum approximately corresponds to E(k) oc k~ 3 . 

This two-component spectrum is a classic feature of 2D turbulence, and is further described 
in reviews of 2D turbulence 



27 



28 29 . Physically the inverse energy cascade corresponds to the 
aggregation of vorticity such that small-scale forcing leads to the growth of large-scale vorticity until 
flows are of the order of the size of the system. At this scale dissipation of some form occurs, possibly 
due to frictional damping at the container's walls. 

Prior to the development of 2D turbulence theory in the 1960s 24 25 [26], Onsager [40{ 



41 



approached turbulent flow in 2D inviscid fluids by focusing on the dynamics of point-like 
centers of vorticity, or point vortices. A collection of point vortices may be constructed to ap- 
proximately represent classical continuum flows, or may quite accurately describe vortex distribu- 
tions in superfluids. In fact, in the same conference proceeding in which Onsager described his 
theory of vortex states in 2D turbulent flows, he also first postulated the existence of quantized 
vorticity in superfluids, and implied that superfluids could be ideal testing grounds for his theo- 
retical arguments of 2D turbulence [lo]. Onsager argued that a system with a large number of 
vortices would maximize entropy when vortices of the same sense of circulation form aggregates, 
giving more room in a finite-volume phase space for the remainder of the vortex cores to dis- 
tribute themselves and move with fewer constraints. Onsager's predictions, which were the first 
hints of the existence of an inverse energy cascade, have yet to be directly confirmed experimentally 
in a superfluid. However, numerical simulations of 2DQT involving forcing in BECs are indeed 
finding evidence for the existence of vortex aggregates that resemble Onsager's predictions [22] 

m. 



2.2. Vortices and 2D Quantum Turbulence in BECs 

Based on the above descriptions of 2D turbulence, direct experimental challenges for the study 
of 2D turbulence in quantum fluids include: (i) reaching a quasi-steady state of forced 2DQT; (ii) 
observation of the aggregation of like-circulation quantized vortices in a superfluid in a state of forced 
2DQT; and (iii) direct measurement of the kinetic energy spectra associated with such turbulent 
states. These challenges all depend on the existence of means to generate 2DQT, so we now turn 
to the subject of vortices in highly oblate BECs, and then focus on describing methods to nucleate 
large disordered distributions of vortices in these compressible quantum fluids. 
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A dilute-gas BEC is often adequately described in terms of an order parameter, or macroscopic 
wavefunction, of the form [l] 

i>{r,t) = ^fj)e i ^ t \ (1) 

where n(r, t) is the atomic density distribution of the BEC and 0(r, t) is the quantum phase profile. 
The dynamics of the wavefunction ip are governed in the mean-field approximation by the Gross- 
Pitaevski equation (GPE) 

Bib h 2 

where g = 47rh 2 a/m characterizes the interaction strength between atoms, a is the s-wave scattering 
length of interatomic interactions, m is the atomic mass, and V (r, t) is the potential due to the trap 
and other external perturbations. 

From a hydrodynamic perspective, a dilute-gas BEC can behave as a tiny droplet of superfluid, 
where the local velocity of superfluid flow is proportional to the local gradient of the phase: 

v = -V(j). (3) 
m 

Consequently, the velocity field of a BEC is irrotational, V x v = 0, except at singularities of the 
quantum phase </>. These singularities are the superfluid-free cores of the vortices, about which fluid 
circulation is quantized 111. A singly quantized vortex is associated with a 27r-phase winding about 
the core, and vortex positions and circulations can therefore be experimentally determined using 



atom interferometry techniques 42 43 44 , or by sequential images that reveal the dynamics of the 
vortices 45 46 47 48 . Most single-shot images of BECs, however, only reveal the location of the 
vortices and not their circulation. 

Vortex cores in a BEC have a diameter on the order of the healing length £ fi] , which is typically 
sub-micron and much less than the wavelength of imaging light. Release of the BEC from the trap 
followed by ballistic expansion of the atom cloud enables the vortices to expand to diameters large 
enough to be optically resolvable, and this method is the standard approach to observing vortices in 
a BEC. In a 3D trap, vortices may also bend with respect to the imaging axis [49], making them even 
more difficult to observe. In two dimensions however, these Kelvin wave excitations are suppressed 
36 . As an example of the clarity that can be obtained in vortex detection, vortices in a highly 



oblate BEC are shown in Fig. [T] the vortices are the circular dark regions within the BEC fluid, 
shown in grayscale with lightness of shade corresponding to integrated column density after ballistic 
expansion from the trap. 

In addition to isolated vortex cores, vortex dipoles may be found in highly oblate or 2D quantum 
fluids. These structures consist of two vortex cores of opposite circulation in close proximity. Vortex 
aggregates, consisting of clusters of vortices of the same circulation, have also been experimentally 
observed in the form of multi-quantum vortex dipoles [46]. In the context of 2DQT in a BEC, 
however, to date only numerical simulations have observed vortex aggregates and their dynamics 
22 . In BECs with a vector order parameter, other types of vortex structures may also be found; 



see Ref. [3] for a summary of vortex experiments in BECs from their initial study in 1999 through 
2010, including vortices in degenerate Fermi gases and multi-component BECs. 

In the hydrodynamic regime, where quantum pressure can be neglected, and with vortices that 
are spaced far enough apart from each other that their core shapes are only very weakly distorted, 
the GPE formalism can be formally linked to the Navier-Stokes equation that plays a key role in 
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Fig. 1. A representative absorption image of an expanded oblate BEC with a resolvable disordered distribution of 
vortex cores. The column density of the expanded cloud is proportional to grayscale, with lighter shades representing 
higher densities. Vortices appear as the dark circular holes within the cloud. 



classical fluid dynamics [38] . By incorporating into the GPE damping due to atomic scattering with 
thermal atoms, a Reynold's number can also be defined [31], although this has yet to be tested 
experimentally; see Ref. [3l] for a discussion of the GPE - Navier-Stokes equation correspondence. 
Given this correspondence, it should not be surprising that the concepts of turbulence may be 
investigated using BECs and the variants of the GPE. 

Much like the vortex tangle that characterizes quantum turbulence in three dimensions [8], we 
would like to similarly physically characterize a state of 2DQT in a BEC. We consider the simple 
picture of a disordered distribution of numerous vortex cores. Such a distribution would satisfy 
the phenomenological criteria of turbulence stated earlier. First, a large disordered distribution of 



vortices would have non-deterministic and chaotic flow dynamics 50 . Second, numerical simulations 
have shown that such distributions can rapidly distribute vorticity 22 3l] or particles (5lJ, even 
if these are initially concentrated in a small region, thereby satisfying the rapid mixing criterion 
of turbulence. Finally, such a distribution of vortices can display a large range of inter- vortex 
separations with no characteristic length scale; in contrast, a vortex lattice has a well-defined length 
scale of the spacing of the vortices in the lattice 31 . We therefore interpret a disordered 2D 
distribution of numerous vortex cores as a phenomenological picture that corresponds to 2DQT. 
Such a picture can be made more quantitative, for example by invoking the methodology and 
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33 



results of Onsager [40] or vortex position statistics, as has been investigated numerically 
[52] . Nevertheless, our general experimental goal is to create states that consist of many vortices 
distributed throughout highly oblate BECs with no discernible or well-defined arrangement. 

As noted above, a primary factor distinguishing 2D and 3D turbulence is the suppression of 
significant flow dynamics in one spatial direction, relative to the flow dynamics in the plane normal 
to that direction. In the case of quantum turbulence, this corresponds to the suppression of Kelvin 
It may then be possible to study 2DQT in 3D BECs, and it should be noted that 
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quasi-2D BECs 53 or even highly oblate BECs are not strictly needed for studies of 2DQT. On the 
other hand, even in highly oblate BECs, vorticity can decay due to vortex-antivortex annihilation. 
It is therefore possible that enstrophy need not be conserved and that a direct energy cascade may 
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appear, as has been recently noted in numerical simulations of decaying 2DQT [30]. In a practical 
sense however, we expect that 2DQT is more readily achievable in highly oblate geometries. Our 
experimental focus is therefore on generating large disordered distributions of vortices in highly 
oblate BECs. 

3. Experimental Setup 




/ \ 



5x imaging system x y 

and CCD camera 

Fig. 2. Experimental configuration for the combined magnetic and optical harmonic trap. A cylindrical lens is used 
to focus a collimated beam of red-detuned laser light into a sheet with a vertical 1/e 2 radius of 23 jum at the location 
of the BEC. The coils for the DC component of the TOP trap, represented as dark gray bars, are located above and 
below the BEC cell. The coils for the AC component, represented as light gray bars, are located on either side of 
the cell normal to the x-axis, and the y-axis (not shown). A lens below the cell collects light for the vertical imaging 
system. 

To obtain the highly oblate 87 Rb BECs desired for our studies of vortex and 2DQT generation 
methods, we use a conventional magnetic time-averaged orbiting potential (TOP) trap [EI] overlaid 
with a 1090-nm red-detuned laser light sheet propagating along the x-axis as illustrated in Fig. [2] 
and Fig. [3ja). A cylindrical lens focuses the laser light in the z direction only, with a waist at 
the position of the BEC such that the laser has horizontal and vertical 1/e 2 radii of (wo y ,wo z ) = 
(2000 /im, 23 /mi). This red-detuned optical potential provides the additional axial (z) confinement 
needed to flatten the BEC while having a minimal effect on the radial (r) confinement. 

We begin our experiments with evaporative cooling of atoms in the \F = 1, mp = — 1) hyperfine 
state of 87 Rb in the purely magnetic TOP trap. Prior to the final evaporation stage, we ramp on 
the trapping laser to a power of ~ 1 W over 4 s, and then cool the system below the BEC transition 
temperature T c ~ 100 nK in the combined magnetic and optical potential. The potential well that 
confines the BEC is an axially symmetric harmonic trap with radial and axial trap frequencies of 
(u; r , uj z ) - 2tt x (8,90) Hz. 

The BECs formed in the combined magnetic and optical potential have up to iV ~ 2 x 10 6 
atoms with a system temperature of T ~ 50 nK and Thomas-Fermi radii of (R r ,R z ) ~ (52,5) jam. 
The chemical potential is /io ~ 8fkj z , and our BECs are consequently far from the quasi-2D limit 
and well into the 3D regime. Nevertheless, the R r : R z ~ H : 1 aspect ratio generated by the 
tight axial confinement suppresses vortex bending and tilting, resulting in a system characterized by 
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Fig. 3. BEC imaging, (a) illustration of the imaging layout. Both horizontal and vertical imaging systems have 
magnification M = 5. Side views of the BEC are obtained with phase contrast imaging along the horizontal imaging 
axis and top views are obtained with absorption imaging along the vertical, (b) a vertically directed absorption 
image of the BEC in the combined magnetic and optical trap immediately after removal of the trapping potential 
but prior to any expansion, (c) a corresponding in situ horizontal phase-contrast image of the BEC in the combined 
magnetic and optical trap, (d) in situ horizontal phase-contrast image of the BEC in the purely magnetic TOP trap 
for comparison, with an aspect ratio of Thomas-Fermi radii of R r : R z = 2 : 1. 



approximately 2D fluid dynamics [36] . Figures [3^b) and [SJc) show images of a condensate held in 
the combined magnetic and optical trap. For comparison, Fig. [3|d) is an in situ image of the BEC 
in the purely magnetic TOP trap with an aspect ratio of R r : R z = 2 : 1. Side views of the BEC 
are obtained with phase contrast imaging along the horizontal (y) imaging axis and views along the 
vertical (z) direction are typically obtained with absorption imaging. 

Vortex cores in our highly oblate BECs have a diameter on the order of the healing length, 
£ ~ 0.4 /am for our system. This is well below the resolution of our imaging system and it is 
currently beyond our ability to detect vortices with in situ images. In order to optically resolve 
the individual vortex cores of highly oblate BECs, we first turn off the TOP trap and allow the 
BEC to expand for 10 ms in the optical potential. The optical trapping potential is then removed, 
and the BEC freely expands for an additional 40 ms. During the entire expansion procedure, a 
separate magnetic field is applied with a gradient along the z direction that cancels the downward 
gravitational force on the cloud. Figure [I] gives a representative absorption image after expansion, 
demonstrating that vortex cores are indeed resolved. In the remainder of this article, all images that 
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show vortices are absorption images taken with the vertical imaging system after a total of 50 ms 
of expansion. 



Optical 
Fiber 



Imagin, 
beam 




TOP trap 
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Fig. 4. Optical layout for the focused blue-detuned Gaussian beam directed axially through the BEC. An image of 
the laser beam can be acquired using the vertical imaging system; an example image is shown. The y push coil (and 
a similar coil along x) can be used to manipulate the position of the center of the magnetic trap. 



In many of the experiments described in the following sections, we make use of an additional, 
focused blue-detuned Gaussian laser beam with a variable power and 1/e 2 beam radius, directed 
axially through the BEC as illustrated in Figure [4] This beam can be used to perturb and excite 
condensates through a variety of methods, which are described below. 

We now turn to various experiments and tests of experimental methods, including relevant re- 
views of previous work in our lab on driving vortices into highly oblate BECs, and discussion of 
previously unpublished results on generation of 2DQT in highly oblate BECs. 



4. Nucleation of Vortex Dipoles and Vortex Clusters 



As described in Neely et al 46 , our group observed the deterministic formation and dynamics of 
vortex dipoles in a BEC by slightly stirring the BEC with a laser beam. In this experiment, vortices 
were nucleated as the BEC was pushed around and past a blue-detuned laser beam that served 
as an impenetrable obstacle. In the frame of the moving BEC, the beam was swept through the 
middle of the BEC, and the vortex dipoles formed in the wake of the laser beam. A vortex dipole 
consists of a vortex and an antivortex - two vortices of opposite circulation - and tends to exist and 
move as a single excitation with linear momentum as long as the two vortices remain far from other 
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Fig. 5. The BEC initial state and experimental sequence, (a) Side-view phase-contrast image and (b) axial absorption 
image of a BEC in the highly oblate trap in the absence of the obstacle. Lighter shades indicate higher column densities 
integrated along the line of sight. Our axial and radial trapping frequencies are measured to be uo z = 2tt x 90 Hz and 
uj r = 2tt x 8 Hz, respectively, (c) BEC initial state with the obstacle located at x s = —20 ~ fim relative to the BEC 
center. (d),(e) The maximum repulsive potential energy of the obstacle is Uo « l.2fio (where /xo ~ 8huj z is the BEC 
chemical potential) and is held constant during evaporative cooling. It is ramped down linearly as the trap translates; 
relative to the trap center, the beam moves from position x s = — 20 ~ /im to x s = 14 ~ jum over a variable sweep 
time t s . The BEC is then held in the harmonic trap for a variable time th prior to expansion and imaging. Figure 
and caption from Ref. [46], T. W. Neely, E. C. Samson, A. S. Bradley, M. J. Davis, and B. P. Anderson, Phys. Rev. 
Lett. 104, 160401 (2010). Copyright 2010 by the American Physical Society. 



vortices or the boundaries of the superfluid. A vortex dipole can be generated from sound, hence a 
moving obstacle can excite a dipole in a compressible superfluid. Similarly the vortices of a dipole 



can recombine, annihilating one another and generating a pulse of acoustic energy; see Ref. 31 for 
a recent discussion of the role and consequences of compressibility in 2DQT. 



4.1. Overview of experiment 

In this experiment, the relative motion between the BEC and the blue-detuned laser beam occurred 
with a constant velocity, as illustrated in Fig. [5] BECs were created in a highly oblate harmonic trap 
as described above in Sec. [3] and shown in Figs. |5ja) and[5jb). The focused blue-detuned Gaussian 
laser beam with a 1/e 2 radius of 10 fim was directed axially through the BEC 20 /im to the left of 
the BEC center as shown in Fig.[5|c). A magnetic bias field was used to translate the harmonic trap 
minimum horizontally along the x-axis at a constant velocity until the laser beam was located 14 
/am to the right of the center of the BEC. At the same time the intensity of the beam was ramped 
linearly from a maximum of Uo ~ 1-2/io down to Uq = as shown in Figs. |5jd) and|5je). Here 
/io ~ Shw z is the chemical potential of the BEC in the highly oblate harmonic trap. At the end of 
the sweep, the BEC was held in the harmonic trap for varying hold times before expanding and 
imaging. 

Neely et al found a critical velocity necessary for observing vortex dipole nucleation of 170-190 
jLtm/s. The repeatability and coherence of this process allowed the dynamics of the vortex dipole to be 
mapped using images of BECs from multiple experimental runs, with each image taken at increasing 
values of t^. As shown in Fig. [6| the vortex dipole orbited the BEC with a period of ~ 1.2 s. GPE 
simulations of the procedure provided visual and quantitative comparisons with the experimentally 
observed vortex dynamics; see Fig.^c). Hall' s group at Amherst College has also observed a variety 
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Fig. 6. Images showing the first orbit of vortex dipole dynamics, (a) Back-to-back expansion images from the 
experiment with 200 ms of successive hold time between the 180-jum-square images, (b) 62-/im-square images from 
numerical calculations of the GPE obtained for conditions similar to the data of sequence (a), but for a temperature 
of T = (i.e., no damping was used). The orbital period is ~ 1.2 s, and the apparent vortex core size is smaller in the 
simulations because we show in-trap numerical data, (c) Black and dark gray small circles show average positions of 
the vortices x v and y v (as fractions of the BEC radius) of each of the two vortices from 5 sequences of experimental 
data identical to that of sequence (a). The larger circle around each average position point represents the standard 
deviation of the vortex positions at that specific hold time, and is calculated from the 5 images obtained at that 
time step. A continuous dipole trajectory from sequence (b) is rescaled to the Thomas-Fermi radius of the expanded 
experimental images, and superimposed as solid lines on the experimental data; no further adjustments are made for 
this comparison. Figure and caption adapted from Ref. [46], T. W. Neely, E. C. Samson, A. S. Bradley, M. J. Davis, 
and B. P. Anderson, Phys. Rev. Lett. 104, 160401 (2010). Copyright 2010 by the American Physical Society. 



of vortex dipole dynamics, using vortices spontaneously trapped in a BEC during the phase transition 



48 



55 and a detection method that permits multiple imaging of vortex positions in a single BEC 47; 



Sweeping the BEC past the barrier at velocities much faster than the first critical velocity gen- 
erated multiply quantized vortex dipoles. At a coarse scale, multiply charged dipoles behave like 
pairs of highly charged vortices of opposite circulation. The loci of vorticity orbits the condensate 
just as in the case of a singly charged dipole, albeit with a shorter period indicative of faster fluid 
flow. A doubly charged vortex dipole, for example, orbited with a period of ~ 0.8 s. As shown in 
Fig. [7| at finer scales these loci of vorticity are actually aggregates of singly quantized vortices of the 
same circulation. Both singly and multiply charged dipoles exhibited lifetimes up to many seconds 
indicating that vortex- ant ivortex recombination may be suppressed and that vortex dipoles may be 
meta-stable in highly oblate BECs. 




Fig. 7. An experimental image of a triply charged vortex dipole in a ballistically expanded BEC. Adapted from 
Ref. [46]. 
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4.2. Implications for the experimental study of 2DQT 

The generation of acoustic energy and the nucleation of vortex dipoles due to a moving obstacle 
are associated with injection of kinetic energy (both compressible and incompressible) into the 
condensate. The nucleation of vortices is also associated with the injection of enstrophy into the 



fluid, with the net enstrophy being proportional to the number of vortices |30||31]. We may view this 
basic laser sweep mechanism as a possible tool for the generation of 2DQT in a BEC, however in the 
experiment described above, the deterministic dynamics of the vortex dipoles indicate that a state 
of 2DQT was not observed. The length scale associated with incompressible energy and enstrophy 
injection is most readily interpreted as the spatial separation between the vortices of the dipole at 
the instant of dipole nucleation; see Ref. [3l] for further discussion of this point. With continuous 
laser stirring within the BEC, one may envision the generation of a 2DQT state, and indeed such 
a state has been observed as described in Sec. ^\ See also the relevant simulations of Refs. [7] and 
[32] that involve an obstacle repeatedly sweeping through a BEC. 

Whether or not a BEC with many vortices and antivortices would generally show traits asso- 
ciated with an inverse energy cascade or a direct energy cascade may strongly depend on specific 
experimental circumstances. In particular, vort ex- ant i vortex recombination is one mechanism for 
dissipation of enstrophy in 2D fluid flow and might inhibit the formation of an inverse energy cascade, 
especially for decaying 2DQT [30]. The balance between vortex generation and vortex annihilation 
may influence observed characteristics of turbulence, such as energy spectra, kinetic energy flux, and 
aggregation of vorticity. 

Finally, the vortex aggregates observed in this experiment are directly generated by the relative 
motion between the laser obstacle and the BEC, and we do not interpret their deterministic formation 



as a consequence of 2D turbulence as proposed by Onsager 40 . Nevertheless, their long lifetimes 
(at least one orbital period) and steady alignment with the tight (z) trapping axis may be taken as 
an indication that vortex aggregates in a highly oblate BEC can be long-lived, and vortex-antivortex 
annihilation rates may be slow enough to reach regimes of constant or increasing net enstrophy. In 
this sense, the observation of vortex aggregates in the experiment of Neely et al gives hope for future 
experimental observation of vortex aggregates in 2DQT BEC experiments. 



5. Generating Turbulent States by Modulating the Magnetic Trapping Potential 

Our TOP trap combines a DC quadrupole magnetic field having variable axial gradient B' z with an 
AC bias field of magnitude Bo that rotates in the x-y plane with frequency ujtop = 27r x 4 kHz. 
The time-dependent bias field with x and y components 



B x = B cos(uj TOP t) (4) 
B y = B sin(cj TO pt) (5) 



results in the minimum of the quadrupole field orbiting in a circle within the x-y plane. For an 
atom with magnetic dipole moment \i and mass m, the time average of the combined magnetic fields 
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results in a harmonic potential with trap frequencies determined by 

1/2 



8mB 



(i - iff 2 



1/2 



(6) 
(7) 



where the factors involving r] = mg/fiB^ account for the vertical sag due to gravity (with acceleration 
g) of the TOP trap minimum from that of a static quadrupole field. Such traps are now in common 
use; see Ref. [56] for a detailed description of TOP traps. Modulations of the magnetic bias field 
and hence the shape of the trapping potential can excite various modes in the BEC. For example, 
by deforming the trapping potential and rotating the spatial deformation axis, surface waves can 
be excited that decay to numerous vortices 12 15 57 . In this section, we describe the vortex 



distributions and their lifetimes that result when highly oblate BECs are subjected to modulations 
of the magnetic component of the trapping field. 



5.1. Symmetric modulation: harmonic trap 

In this experiment, we modulated the amplitude of the magnetic bias field Bq in Eqs. Q and (J5| 
by adding a small sinusoidal modulation to Bq of the form 

B (t) = B + B mod sin(uj mod t). (8) 

Here B mod is the amplitude of the bias field modulation and cj mo d is the frequency of the modulation. 
This results in a small sinusoidal oscillation of the radial trap frequency u Y about the static value, 
but the trap remains axially symmetric. The axial trapping frequency u z is primarily determined by 
the optical potential and therefore remains approximately constant. Our aim in this experiment was 
to determine if a weak symmetric modulation would drive vorticity into the BEC such that there 
was approximately no net angular momentum transfer, i.e., such that if any vortices were generated, 
approximately equal numbers of vortices of both circulations would be created. 

In our procedure, we began with a highly oblate BEC in the combined magnetic and optical 
harmonic trap. We then modulated the amplitude of the rotating magnetic bias field at frequency 
w mo d = 2tt x 9 Hz and amplitude B mo( ± ~ 0.05 Bo for varying modulation time £ mo d- After the 
modulation we held the BEC in the static harmonic trap for varying hold times th, then expanded 
and imaged, looking for the presence of vortices. 

Figure [8ja) shows images taken at increasing modulation times from t mo a = 0.2 s to 1.5 s. 
Vortex cores started to appear at t mo( ± = 0.4 s and by t mo d = 0.6 s we observed disordered vortex 
distributions. Vortex cores appear to be generated at the outer edge of the condensate and make 
their way into the center. In the next step of the experiment, we modulated the trap as described for 
t mo d= 1.5 s before removing the modulation and returning the BEC to the unmodulated harmonic 
trap. We added increasing post-modulation hold times th and then expanded and imaged the BEC. 
As shown in Fig. [8^b) the vortex distribution became better resolved with increasing th and decayed 
to states with only a few cores remaining after th = 2 s. By th = 8 s no vortex cores remained. The 
condensate was completely destroyed for £ mo d = 10 s. 
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Fig. 8. Vortex distributions after varying modulation and hold times, (a) images taken at increasing modulation 
times from £ mo d — 0-2 s to 1.5 s. Vortex cores start to appear at t mo d — 0-4 s and and by t mo d — 0-6 s we see a 
disordered vortex distribution, (b) images taken for increasing hold times after modulating for £ mo d — 1-5 s. 



5.2. Symmetric modulation: annular trap 

In a similar experiment, we modulated the trapping potential of a highly oblate annular trap. To 
do this, we modulated the strength of the rotating bias field of the TOP trap in the same manner as 
described in Sec. |5.1| We added to the harmonic trap a focused blue-detuned Gaussian laser beam, 
directed axially through the BEC in the configuration shown in Fig. |4j Example images of such a 
BEC are shown in Figs. [9ja) and|9jb). Prior to imaging, the blue-detuned laser barrier was ramped 
off slowly enough to avoid significantly perturbing the condensate, allowing the BEC to expand 
from a harmonic trap for imaging. With cj mo d = 2tt x 6.5 Hz and £? m od ~ O.LBo, numerous vortices 
were quickly nucleated over just a few ms, a significant change from the harmonic trap. As shown 
in Fig. [9jc) vortex cores were observed to be distributed throughout the BEC, with the first cores 
appearing at least an order of magnitude earlier than in the purely harmonic trap. Compared to 



the results of Sec. 5.1 , many more cores were visible using this modulation technique in the annular 
BEC, and the cores have higher contrast. We speculate that the presence of the blue-detuned barrier 
acted as a vortex dipole nucleation site located within the BEC, and that the vortex dipoles quickly 
dissociated and formed a state of 2DQT. Again, the modulation process was designed to be axially 
symmetric apart from experimental errors in beam positioning and trap roundness, and should have 
nominally added no net angular momentum to the condensate. 

With additional modulation time, the BEC was driven into a state of high excitation in which 
individual vortex cores were not easily resolvable or countable. Such a state is shown in the rightmost 
image ofFig.gc). This process shows that too much modulation drives the BEC into a highly excited 
far-from-equilibrium state, a signature that it may be possible to identify a weaker modulation 
strength such that vortex nucleation rates and vortex decay rates due to annihilation and thermal 
damping at the system boundary would be balanced for a given temperature. Finding this optimum 
balance is a goal of future research. 
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Fig. 9. Modulating the radial trapping frequency for a BEC in an annular trapping geometry, (a) side view of 
an annular BEC imaged using in situ phase contrast imaging, (b) top-down view of an annular BEC imaged using 
absorption imaging immediately after turning off the trapping potential. For these two images, the blue-detuned beam 
provides a barrier of height well above the chemical potential, confining the atoms to an annular trap. The annular 
trap used for the sequence of images in (c) had a narrower, weaker central barrier, (c) images taken at increasing 
modulation times with no hold after modulation. The first image shows the expanded BEC for t mo d = s to show 
that the hole created by the blue-detuned beam has completely filled between beam ramp down and imaging. For the 
annular trapping geometry vortices appear after much shorter modulation times than for the harmonic trapping case. 



5.3. Rotation of an elliptical magnetic trapping potential 

The straightforward experimental method of expansion imaging used in our experiments does not 
easily permit measurement of the circulation of vortices. While the interferometric or dynamic 
methods mentioned above can be used for this purpose, these have not yet been applied to turbulent 
states of BECs. In order to look for the clustering of vortex cores of like-circulation, an alternative is 
to examine statistics of vortex distributions containing only large numbers of cores of identical circu- 
lation. One possible approach to generate such distributions is by rotating the trapping potential, as 



has been utilized in numerous experiments of BECs with different aspect ratios, such as Refs. 12 15 



57 ; see also Ref. [3] for an overview of numerous other experiments that utilized this technique. Our 
approach follows such previous work, and is implemented by squeezing and rotating the magnetic 
trapping field as first implemented by Hodby et a/ [57]. 

In order to rotate and add angular momentum to the atomic cloud it is necessary to break the 
symmetry of the trapping potential. In our case this is done by adding an additional bias field to 
the one described by Eqs. Q and (J5| so that now 

£ x = B cos(cjtop^) + B e cos(u mod t) (9) 
B y = B sm(u TOP t) - B e sm(u mod t) (10) 

where B e ~ 0.1 Bo, and u; mo d is the frequency of the modulating bias field. If u mo a = cjtop the 
bias field causes the time-averaged trap to have stationary ellipsoidal potential energy surfaces in 
the x-y plane with the ratio of the minor and major axes determined by B e and Bq; f° r °ur case 
this ratio is ~ 0.8. Increasing B e with respect to Bq increases this ellipticity. If oj mo d ^ oj^op the 
ellipse rotates in the x-y plane with frequency uj s — |co> mo d — cjtop|/2. I n the TOP trap, choosing 
uj s ~ 0.7co> r excites the quadrupole mode which decays to a collection of quantized vortex cores of 
the same circulation. The minimum energy configuration of the rotating BEC consists of a large 
regular lattice of vortices [12]. A representative vortex lattice formed by spinning the BEC in our 



purely magnetic TOP trap (i.e., without the red-detuned trapping potential) is shown in Fig. 10 

Our investigations of highly oblate BECs subjected to trap deformation and rotation started with 
a BEC in the combined magnetic and optical harmonic trap. We then applied the magnetic field 
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Fig. 10. An image of a slightly irregular vortex lattice in a rotating BEC released from our TOP trap. 

ellipticity and spun in the highly oblate harmonic trap for time t s at frequency uj s with ellipticity 
B e . After spinning we returned to the symmetric harmonic trap and held for time th while the BEC 
shape deformations damped out. Finally, we let the condensate expand and we imaged the cloud. 




Fig. 11. Comparison of number of vortex cores versus spin time for different ellipticities. All data points were 
generated by spinning at uj s = 2tt x 6 Hz for a variable spin time t s and then holding for 2.5 s after the spin prior to 
expansion. Open circles represent a single set of data taken with B e = 0.05 £?o, gray squares represent B e = 0.07 Bo, 
and black triangles represent B e = 0.09 So- Both the peak number of vortex cores and t s corresponding to the peak 
number of vortex cores seem to be dependent on ellipticity but all three ellipticities show decreasing vortex number 
for longer spin times. 

By varying spin frequencies, we found a resonance condition for generating large numbers of 
vortices at uj s = 2tt x 6 Hz. Here uj s ~ 0.75co> r is consistent with previous observations of exciting 
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the quadrupole mode in 3D harmonic traps at co> s ~ 0.7cj r [16] . We measured the number of vortex 
cores versus spin time for different ellipticities at this resonance frequency. Figure 11 shows three 
sets of data corresponding to ellipticities of B e — 0.05 £?o, 0.07 Bo and 0.09 Bo. All data points were 
generated by spinning for a variable spin time t s and then holding for th = 2.5 s after the spin and 
before expansion. Both the peak number of vortex cores and t s corresponding to the peak number 
of vortex cores appear to be dependent on ellipticity with peak vortex number occurring at a much 
shorter spin time for B e — 0.09 Bo and 0.07 Bo than for B e = 0.05 Bo. All three ellipticities show 
decreased vortex numbers for longer spin times on the order of t s = 8 s indicating that vortex cores 
may be annihilating or leaving the system before crystalizing into a lattice. 
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Fig. 12. 350-/im-square expansion images of an oblate BEC for varying t s (indicated above images) at increments 
of 100 ms. Spin frequency and ellipticity were held constant at uj s = 2tt x 6 Hz and B e = 0.07 Bo. (a) and (c) BECs 
were expanded and imaged directly after spinning, t h = 0. (b) and (d) BECs were held for 400 ms after spinning, t^ 
= 400 ms, then expanded and imaged. 



Figure 12 visually tracks the evolution of the condensate during spinning. Figures 12 'a) and[l2[c) 
show successive images at 100 ms intervals for t B = 0.1 to 1.0 s. The BECs were imaged directly 
after spinning with no hold in the axially symmetric harmonic trap (B e = 0) prior to expansion. 
Figures 12 'b) and 12 'd) show successive images at 100 ms intervals for t s = 0.1 to 1.0 s with an 
additional 400 ms hold time in the axially symmetric harmonic trap prior to expansion and imaging. 
The hold time in the symmetric harmonic trap seems important for the nucleation of vortices as 
vortex cores appear after t s = 400 ms and th = 400 ms, but do not appear until t s = 1.0 s for the 
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case of th = 0. 
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Fig. 13. 350-/im-square expansion images of vortex distributions in the highly oblate BEC for varying t^. Spin time, 
frequency, and ellipticity were held constant at t s = 1.0 s, w s = 27r x 6 Hz and B e = 0.07£?o- Images were taken for 
variable th from to 2.5 s. Multiple images for = 2.5 s are shown to give a representative sample of the variation 
in vortex distribution. 
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Fig. 14. 350-/xm-square expansion images of vortex distributions in the highly oblate BEC for the varying hold times 
th shown from 3 to 8 s. By this time, BEC shape deformations have largely damped out. Spin time, frequency, and 
ellipticity were held constant at t s = 1.0 s, uu s = 2tt x 6 Hz and B e = 0.07 Bq. 



Figures [T3| and [14] show additional distributions of vortex cores nucleated by spinning and relaxing 
over subsequent hold times. Here t s = 1.0 s is held constant and the BECs are imaged after increasing 
values of t^. Vortices appear to be nucleated on the outer edge of the BEC. Presumably there is net 
angular momentum added to the condensate, so we expect the vortex cores to be predominately of 
the same circulation, at least once the system has relaxed to a meta-stable configuration. However, 
it may be the case that the net circulation is large, but that numerous vortices of the opposite 
circulation are also generated and present in these images, at least prior to the point at which a 
disordered distribution is uniformly distributed throughout the BEC. If vortices are indeed of the 
same circulation in images such as the ones in the bottom row of Fig. [l3j then such states may be 
candidates for experimental measurement of vortex power-law distributions in a BEC, as has been 
analytically described for homogeneous BECs [31]. 
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As observed, large, disordered vortex distributions can be generated by exciting collective modes 
of the BEC. The shape deformations damp more quickly than the number of vortex cores, leaving 
open the possibility of finding parameters where studies of decaying 2DQT could be performed before 
the system spins down to a state with no vortices. Spinning introduces a net angular momentum 
into the condensate and should result predominately in vortex cores of the same circulation, making 
this a possible system in which to observe vortex aggregates without time-resolved dynamics mea- 
surements. Nevertheless, this vortex excitation technique does not appear to satisfy the particular 
goal of continuous forcing, although it appears to be a candidate for studies of decaying 2DQT. It is 
interesting to note that even for long hold times we find no evidence of a vortex lattice in the highly 
oblate BEC. Thermalization times for a lattice in a highly oblate BEC may be beyond the lifetime 
of our BECs; see Ref. 58 for further discussion of this issue. 



6. Generating Turbulent States with a Stationary Blue-Detuned Laser 

In the experiments described in this section, we investigated the response of a BEC to time-dependent 
perturbations of the intensity of a focused blue-detuned laser beam that pierced the BEC. In all 
cases in this section, the relative position between the beam and the BEC was stationary. Our aim 
in these experiments was to locally excite the BEC as an empirical probe of the existence of thermal 
counterflow [59} |60] in a region where the BEC was locally depleted of atoms. While vortices were 
observed in all methods examined, the mechanisms for vortex nucleation remain unclear and merit 
further experimental and numerical investigation. 



6.1. Short pulse of blue- detuned laser light 

In this experiment we nucleated vortices by subjecting the BEC to a short pulse of blue-detuned 
laser light. We began by forming a highly oblate BEC in the purely harmonic trap. We then 
instantaneously turned on a blue-detuned Gaussian beam for a pulse time of 7 ms, short compared 
to the radial harmonic oscillator period of 125 ms. The focused blue-detuned Gaussian beam had 
a 1/e 2 radius of 10 fim and was directed axially through the center of the condensate, as shown in 
Fig. [4) The power in the blue-detuned beam was chosen such that the optical potential generated 
by the beam was approximately equal to the chemical potential of the BEC. After the laser pulse 
we held the BEC in the harmonic trap for varying hold times th before expansion and imaging. 



Figure 15 [a) shows an in situ absorption image of the BEC with the blue-detuned beam on, clearly 
penetrating the condensate. 

The laser pulse generated a shock wave that can be observed as the central high density region 
in the first image in Fig. [T5^b). Subsequent images show the condensate for longer hold times. 
Vortices appear to enter the condensate from the outer boundary as the shock wave propagates. 
Vortices eventually leave the system, most likely through vort ex- ant i vortex recombination and ther- 
mal damping. Apart from the shock wave, this perturbation created vortices with little residual 



excitation of the BEC. Hoefer et al 61 performed experimental and numerical studies of the dis- 
persive Shockwaves generated when a BEC was subjected to a short pulse of laser light, however 
their experiments were performed with a prolate BEC and they do not report observations of vortex 
cores generated as a result of the laser pulse. 



20 



K. E. Wilson, E. C. Samson, Z. L. Newman, T. W. 



Neely, and B. P. Anderson 



(b) 



(a) 


















■ l 


1 f 


EJ 








i i 




50|im 






1 - 1 




t h = 10 ms 


t h = 25 ms 


t h = 50 ms 


t h = 125 ms 



t h = 250 ms 



t = 500 ms 



t*= 1 s 



t, = 2s 



Fig. 15. Vortices generated by subjecting the BEC to a 7 ms pulse of blue-detuned laser light directed axially through 
the center of the condensate, (a) in situ absorption image of the BEC with the blue-detuned beam penetrating the 
center of the BEC. The beam intensity is chosen such that the beam potential is approximately equal to the condensate 
chemical potential, (b) Images show successive hold times after the 7ms laser pulse. 



6.2. Intensity modulation of a blue-detuned laser beam 

In this experiment we formed a highly oblate BEC in the purely harmonic trap, then turned on and 
sinusoidally modulated the intensity of a focused blue-detuned Gaussian beam for a varying time 
^mod- After the modulation we held the BEC in the purely harmonic trap for varying time t^, then 
expanded and imaged. The blue-detuned beam had a 1/e 2 radius of 10 jam and was directed axially 
through the center of the condensate, as shown in Fig. [4] During the modulation time the optical 
potential U (t) generated by the beam followed 

U(t) = U sin 2 (owit/2) (11) 

for < t < tmocb with maximum repulsive potential energy Uq, and frequency cj mo d. For all other 
times, the beam was turned completely off. We fixed the modulation time at integer multiples of 
the modulation period, r mo d = 27r/u; mo d, such that the intensity of the blue-detuned beam always 
started and ended at zero, and we did not have to be concerned with ramping off the blue-detuned 
beam for imaging. As in the laser pulse technique, the blue-detuned beam acted as a perturbation 



to the confining potential. Figure [16] shows a sequence of images for varying values of after 
modulating for t mo( ± = 187.5 ms, with Uo ~ 0.3/io (where uo ~ 8So; z ), and frequency u mo( ± = 2tt x 16 
Hz. This was approximately twice the radial trap frequency. An in situ image of the BEC in the 
harmonic trap with the blue-detuned beam aligned in the center is shown in the leftmost image of 



Fig. 17 



While this method generated a large number of vortex cores, we also observed large-scale breath- 
ing oscillations in the radial dimension of the condensate. The radius of the expanded cloud oscillated 
between ~ 40 urn and ~ 180 urn with a period of ~ 70 ms and eventually damped out with an expo- 



nential decay time constant Td am p °° 500 ms. As shown in Fig. 17, a disordered distribution of vortex 
cores remained in the condensate even after the bulk oscillations subsided, with ~ 8 cores remaining 
for th = 2 s. It is not entirely surprising that we induced bulk fluid oscillations, given that we were 
forcing the BEC at a frequency that was twice the radial trapping frequency. However, further study 
is needed for positively identifying the vortex generation mechanism in this experiment. 

In the sequence of images shown in Fig. [16} it is unclear where vortex nucleation occurs. To 
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Fig. 16. 250-/im-square absorption images acquired after the hold times th indicated. The condensate undergoes 
large breathing oscillations in the trap, and this oscillation leads to a periodic variation in expanded BEC radius 
between ~ 40 jim and ~ 180 jim. 
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Fig. 17. 250-/im-square absorption images acquired after the hold times indicated. The first image is an in situ 
image of the BEC in the harmonic trap with the blue-detuned beam partially penetrating the condensate; note that 
the optical potential strength is well below the BEC chemical potential. 



further explore the nucleation mechanism, we aligned the blue-detuned beam near the edge of the 
BEC, and modulated for t mo( ± = 62.5 ms, equivalent to one sinusoidal pulse of the laser light, with 
Uo ~ 0.3/io, and frequency u; mo d = 27r x 16 Hz. After modulation, we held the BEC in the purely 
harmonic trap for th = 40 ms, then expanded and imaged the BEC. As shown in Fig. [l8j the vortex 
cores that resulted from the modulation appear to form near the location of the focused laser beam. 

Lastly, we replaced the axially symmetric (wq x = wq y ) focused blue-detuned Gaussian beam 
with a blue-detuned light sheet, focused along the x-axis with a 1/e 2 radius wq x = 10 /im, spatially 
extended along the y direction, and directed axially (along z) through the center of the BEC. The 
width of the beam along the ?/-axis was much larger than the diameter of the condensate so that 



beam extended beyond the edge of the condensate in the ^/-direction. Figure 19 'a) is a vertical 



absorption image of the unexpanded highly oblate BEC with the elongated blue-detuned beam 
partially penetrating the condensate. We modulated at oj mo d = 27r x 16 Hz for t mo d = 62.5 ms and 
Uo ~ 0.3/io, held for th, then expanded and imaged the BEC. Figures 19 'b) and [l9^c) were taken 
after th = 40 ms. Here the vortices appear to be nucleated along the long axis of the elongated 
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Fig. 18. Localization of vortex cores at the location of the blue-detuned beam. All images are 250-/im-square 
absorption images. The top row of images were taken with the BEC in the harmonic trap and the maximum height 
of the blue-detuned potential at Uo ~ 0.3/io- The locations of the laser beam correspond to the positions designated 
by the arrow. The bottom row of images were acquired after = 40 ms followed by the expansion procedure. Note 
the correlation between beam position and the position of the vortex cores in each vertical pair of images. 



beam. The vortex cores shown in these images are not completely resolved but we suspect that the 
vortices are being nucleated as dipoles in a similar manner to the breakdown of a soliton due to the 
snake instability in a BEC [62 63 64 . Alternatively, these features may be acoustic precursors to 
vortex dipole formation [65] . 

As with the axially symmetric blue-detuned potential, we observed bulk excitations in the con- 
densate. Figure |20| shows a sequence of images taken with cj mo d = 2n x 16 Hz, t mo( ± = 62.5 ms, and 
Uo ~ 0.5/io, for varying hold times t^. The oscillations in x and y are now out of phase by ~ 90°. 
Vortex cores appear to be generated at the location of the beam and move to the outer boundary. 
In particular the images corresponding to < 40 ms show vortices aligned with the long axis of 
the beam, but by = 50 ms the vortices are located along the outer boundary of the condensate 



(a) (b) (c) 

Fig. 19. Localization of vortex cores at the location of the blue-detuned beam. All images are 250-/xm-square 
absorption images, (a) taken with the BEC in the harmonic trap and the blue-detuned potential on at a strength of 
Uo ~ 0.5/io- in (b) and (c), expansion images were taken after = 40 ms. For these images Uo ~ 0.3/io- Note that 
the vortex cores are localized along the long axis of the blue-detuned beam and that they appear to be nucleating in 
pairs. 
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Fig. 20. 250-/im-square absorption images acquired after the hold times t^ indicated. The condensate undergoes 
bulk oscillations that are out of phase in x and y. Vortex cores appear to nucleate at the location of the beam and 
then move throughout the BEC. 

with no cores in the center. Later images show more disordered distributions of cores. Again, the 
mechanism for nucleation is not clear. 

In all of these methods, intensity modulations of a blue-detuned laser beam were observed to be 
effective for nucleating vortices in a BEC. Although the mechanisms for vortex nucleation remain 
unclear, we observe correlations between the position of the laser beam and the site of vortex 
generation within the BEC. With careful parameter selection, it may be possible to use this method 
for controlled vortex generation rates. However, one must also be careful not to significantly excite 
shape oscillations, as these would make studies of 2DQT difficult. Perhaps with further adjustment 
of the modulation rates and times, shape oscillations could be minimized, and methods of this sort 
could be utilized for 2DQT studies. 

7. Stirring with a Blue-Detuned Laser Beam 

Rather than using intensity modulation of a blue-detuned laser beam, the experiment described in 
Neely et al [22] used small-scale stirring of the BEC with a blue-detuned Gaussian beam to generate 
disordered vortex distributions in a highly oblate annular trapping potential. It was also observed 
that vortices coalesced into large-scale flow in an annular trap, but here we review only the vortex 
distributions observed via stirring. 

This experiment was performed with the BEC in an annular trap created with a focused blue- 
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Fig. 21. Timing sequence used to study 2DQT. (a) and (b), experimental parameters vs. time. In (a), the radio 
frequency (RF) field is adjusted so that the RF value is jumped further away from cutting into the BEC while the 
BEC is stirred with the laser beam. After the stirring period and a 1.17-ms hold, the RF cuts into the cloud to 
lower the temperature of the system, (b) illustration of the sequence of events during BEC stirring, (c) and (d), 
experimental in situ column-density images of the BEC immediately prior to the stir, as viewed (c) in the plane of 
2D trapping and (d) along the z axis. Lighter grayscale shades indicate larger column densities, (e) illustration of 
stirring, the arrow shows the trajectory of the harmonic trap center relative to the larger fluid-free region created 
by the laser barrier, (f) in situ image of the BEC 10 s after stirring; vortices are not observable, necessitating an 
expansion stage to resolve them. 



detuned Gaussian beam, with a 1/e 2 radius of 23 /xm, directed axially through the center of the 
BEC. As shown in Fig. 21 the blue-detuned beam penetrated the BEC with a barrier height of 
Uo ~ 1.5/io where /xq ~ 8hu z . Prior to stirring, the BEC was held at a temperature T ~ 0.9Xc, 
where Tq ~ 116 nK was the critical temperature for the BEC phase transition. 

At t = magnetic bias coils were used to move the center of the harmonic trap in a 5.7-/xm- 
diameter, off-center circle about the stationary, blue-detuned barrier. At t = 0.333 s, at the end of 
the stirring motion, the center of the harmonic trap again coincided with the blue-detuned beam, 
and the BEC was held in this annular trap for varying hold times up to ~ 50 s. At = 1.17 s 
the temperature of the BEC was reduced to ~ 0.6T C to decrease thermal damping rates and vortex- 
antivortex recombination. At the end of the hold time the blue-detuned beam was ramped off over 
250 ms and the BEC was expanded and imaged. 



Experimental images acquired at varying hold times are shown in Fig. 22 Immediately after 



stirring, a large disordered distribution of vortex cores was observed. By = 0.33 s, vortices had 
begun to coalesce on the central potential barrier, as indicated by the large-vorticity hole in the 
center of the images. Large-scale superflow developed over increasing hold times, as indicated by 
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Fig. 22. 200-/im-square experimental column-density absorption images acquired at the hold times £h indicated above 
the images. Each BEC undergoes ballistic expansion immediately after the central barrier ramp-down in order to 
resolve the vortex cores. Each image is acquired from a separate experimental run. 

the growing large area of vorticity at the center of the BEC, and was observed to persist up to 50 
s. The initial vortex distribution generated by the stirring varied significantly from shot to shot but 
consistently evolved into a large-scale flow. 

The damped GPE was also used to numerically simulate the BEC stirring and subsequent dy- 
namics; see Ref. [22] for details. In addition to the evolution of the large-scale flow observed in both 
experiment and numerics, analysis of the numerical simulations yielded two other results indicative 
of 2D turbulence. First, coherent vortex structures consisting of two cores of the same circulation 
intermittently formed and lasted for long timescales, tantalizing evidence for Onasger-type vortex 
aggregates. Second, analysis of the incompressible kinetic energy spectrum gave E(k) oc /c -5 / 3 for 
k < hi and E(k) oc k~ 3 for k > k{ where k{ is the wavenumber associated with injection of ki- 
netic energy into the BEC. In this case the injection mechanism was energy transfer between the 
sound field associated with the compressible component of the fluid excited by stirring, and the 
vorticity field associated with the incompressible component. Further relationships between stirring 
mechanisms and observed vortex distributions and dynamics are found in Refs. [22] and [31]. With 
the development of new experimental techniques, we hope that such information will soon also be 
experimentally discernible. As it stands, such small-scale stirring with a blue-detuned laser beam 
appears to be a promising mechanism for further studies of 2DQT, perhaps even with multiple beam 
sites or with longer stirring times. 



8. Conclusions 

We have found that it is relatively straightforward to nucleate large, disordered distributions of vortex 
cores in highly oblate BECs, that is, to generate 2DQT in a BEC. We interpret the disordered nature 
of these distributions as analogous to the vortex tangle characteristic of 3D quantum turbulence. 
Many of the vortex distributions that we have observed are long-lived, making them potentially useful 
for studies of 2DQT. However, methods such as modulating a blue-detuned beam in the center of 
the BEC excite bulk modes of the fluid where the entire condensate undergoes shape oscillations. 
In many cases, the vortex nucleation mechanisms are not completely clear. Possible mechanisms for 
vortex generation in the laser modulation experiments include counterflow between the thermal and 



superfluid components [9] [59 60 , vacuum bubble cavitation as described by Berloff and Barenghi 
66 , or perhaps other nonlinear fluid dynamics processes at the boundaries of the quantum fluid. 
While further numerical and experimental investigations are necessary to determine the origins of 
the vortices, we envision many of these highly disordered vortex distributions as starting points for 
further studies of 2DQT. 
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Although we have described a number of vortex generation mechanisms, a primary challenge 
for experimentalists is to characterize 2DQT in a BEC, and the characteristics observed may well 
depend on the vortex generation mechanism used. Experimental observation of an inverse energy 
cascade, vortex aggregation, and kinetic energy spectra are of primary interest. Such characteristics 
might or might not appear via the methods discussed in this article; this remains to be determined. 
Vortex dynamics and circulation measurements will likely prove to be among the most significant 
developments remaining to be developed for 2DQT studies in a BEC. Understanding the temperature 
dependence, and more generally the role of dissipation, are also of high importance. Finally, reaching 
a regime of quasi-steady-state 2DQT is a highly desired goal, one that may well be reached utilizing 
methods described in this article, presumably with a suitable balance of dissipation. 

Although many of the methods we describe in this article represent only the beginnings of investi- 
gations into new methods for the study of 2DQT, we hope that the descriptions of our observations so 
far inspire the development of new experimental and numerical studies of vortex generation mech- 
anisms. Further experimental and theoretical work in the field of 2DQT promises to provide an 
exciting new research direction for compressible quantum fluids. 
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